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Abstract 

A numerical and analytical study of the role of exponentially truncated Levy flights in the 
superdiffusive propagation of fronts in reaction-diffusion systems is presented. The study is based 
on a variation of the Fisher-Kolmogorov equation where the diffusion operator is replaced by a 
A-truncated fractional derivative of order a where 1/A is the characteristic truncation length scale. 
For A = there is no truncation and fronts exhibit exponential acceleration and algebraic decaying 
tails. It is shown that for A 7^ this phenomenology prevails in the intermediate asymptotic regime 
{xt) 1 ^ <C £ ^ 1/A where x is the diffusion constant. Outside the intermediate asymptotic regime, 
i.e. for x > I/A, the tail of the front exhibits the tempered decay 4> ~ e~ Xx / 'x^ 1+a ^ , the acceleration 
is transient, and the front velocity, vl, approaches the terminal speed v* = (7 — A a x)M as * ~~ * °°> 
where it is assumed that 7 > X a x with 7 denoting the growth rate of the reaction kinetics. However, 
the convergence of this process is algebraic, vl ~ v* — a/ (At), which is very slow compared to the 
exponential convergence observed in the diffusive (Gaussian) case. An over-truncated regime in 
which the characteristic truncation length scale is shorter than the length scale of the decay of the 
initial condition, 1/u, is also identified. In this extreme regime, fronts exhibit exponential tails, 
cf) rvj e _ra , and move at the constant velocity, v = (7 — A a x) l v - 
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I. INTRODUCTION 



Reaction-diffusion systems have played a predominant role in the study of pattern for- 
mation and nonlinear dynamics in a large class of phenomena of interest to physics, biology, 
chemistry and engineering, see for example Refs. and references therein. One of the 

simplest reaction-diffusion systems is the extensively studied Fisher-Kolmogorov model that 
describes the dynamics of a scalar field, 0, in a one- dimensional domain, 

^ = x<9^ + 70(1-0) , (1) 

where x denotes the diffusivity and 7 is a constant. The nontrivial dynamics of reaction 
diffusion systems in general stems from the competition between the diffusivity and the 
non-linearity. In the case of the Fisher-Kolmogorov equation this competition leads to the 
propagation of fronts in which the stable, = 1, state advances through the destabilization 
of the = unstable state. 

An important and often overlooked assumption in reaction-diffusion models is the use 
of Laplacian operators, x^ 2 , for the description of transport. The use of these operators 
is motivated by the Fourier-Fick's model according to which the flux q is assumed to be 
proportional to the gradient of the concentration, q = — %V0. This local prescription, 
together with the conservation law, d t (j) = — V q, leads to the Laplacian, diffusive transport 
operator. From the statistical mechanics point of view this prescription is linked to the 
assumption that the underlying "microscopic" transport process is driven by an uncorrelated, 
Markovian, Gaussian process. However, despite its relative success, experimental, numerical 



and theoretical 
example Refs 



evidence indicates that the diffusion model has limited applicability, see for 



• 2! Q, 5), lf| and references therein. Therefore, a problem of considerable interest 
is the study of the role of anomalous diffusion, and Levy flights in particular, in reaction 
diffusion systems. 

y work on reaction-anomalous-diffusion systems include Refs. [lOj. Refer- 



Ear 
ce|? 



ence |7| studied bistable reaction processes and anomalous diffusion caused by Levy flights. 
The interplay of sub-diffusion and Turing instabilities was discussed in Ref. [8j. The role of 
super diffusive transport in the acceleration and algebraic decay of fronts was studied in the 
context of a probabilistic approach in Ref. [9j and in the context of an equivalent fractional 
Fisher-Kolmogorov equation in Ref. [k]] . It is interesting to note that fronts in chaotic couple 
map lattices with long-range couplings exhibit an analogous phenomenology as discussed in 



Ref. 



111 ]. More recent works include: the study of analytic solutions of fractional reaction- 



diffusion 



121 ]; the study of a reaction-diffusion system with a bistable reaction term and 



directional anomalous diffusion 



equations [yj] ; the study of Turing instabilities [lf| ; the study of the effect of superdiffusion 



131 ]; the study of the construction of reaction-sub diffusion 
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; and the study of the fractional 



181 ] among others. 



on pattern formation selection in the Brusselator model 
Ginzburg-Landau and Kuramoto-Sivashinsky equations 

Evidence of Levy flights has been found in laboratory experiments, simple models and 
numerical studies of turbulent transport and the use of fractional diffusion models to describe 
these problems has been well-documented in the literature 5j. However, it is plausible that 
finite-size domain and decorrelation effects (among other effects) might have an impact 
on the Levy flight dynamics. The evaluation of the role of these "non-ideal" effects on 
reaction-anomalous-diffusion systems is a problem of considerable practical relevance. Of 
particular importance is to determine how, and to which degree, these effects might mask the 
underlying Levy statistics. The effect of fluctuations caused by finite number of particles per 
volume on the super diffusive propagation of fronts was studied in Ref. [15]. Here we focus 
on the role of truncation effects on Levy flights driving superdiffusive front propagation. 

Asymptotic analysis plays an important role in the evaluation of non-ideal Levy flight 
effects. In particular, it is quite possible that because of non-ideal effects the statistics of 
the system will eventually converge to Gaussian. However, the key issue is to determine the 
duration of the non-diffusive transient and the rate of convergence to Gaussian statistics. 
This point is clearly illustrated in the ultraslow convergence to Gaussian statistics in the 
presence of truncated Levy flights. In this case it has been observed that although the statis- 
tics eventually converges to Gaussian (because of the central limit theorem) a remarkably 
large number of steps is needed, and therefore the system effectively behaves non-diffusively 



in the intermediate asymptotic regime of practical interest [19|, [20), |21|, |22j, |23J . In the present 
paper we explore to which degree a similar situation occurs in the case of front propagation. 
In particular, in Ref. [lO] it was shown that in the fractional Fisher-Kolmogorov equation 
fronts decay algebraically and exhibit exponential acceleration. The goal of this paper is 
to present a numerical and analytical asymptotic study of the effect of truncation on these 
phenomena. One problem of special interest is to determine if there is an intermediate 
asymptotic regime where the effects of truncation are negligible and where the fronts accel- 



erate and exhibit algebraic tails. Outside such intermediate asymptotic regime it is expected 
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that the truncation effects will become dominant and that as t — > oo the front dynamics will 
eventually approach in some sense the diffusive Fisher-Kolmogorov dynamics. However, the 
key issue is how long will this take. This bring us to the second problem of interest in this 
paper which to determine the rate of convergence to the constant velocity and exponential 
tails characteristic of the diffusive front propagation regime. 

Our approach is based on the use of truncated fractional diffusion operators. Fractional 
derivatives provide a powerful framework to model non-diffusive transport processes 
These operators incorporate long-range, non-local transport through the use of slowly de- 
caying kernels. In particular, in fractional diffusion the Laplacian is replaced by an integral 
operator of the form d 2 x J <f>(x' ,t)K(x — x')dx', where the kernel K has the asymptotic be- 
havior, K ~ x~ a+1 , for 1 < a < 2. In the context of the Continuous Time Random Walk 
model the exponent a is related to the stability index of the underlying Levy process jfj]. 
However, a potential drawback of this description is that a-stable Levy processes have di- 
vergent second moments because their corresponding densities decay as ~ x~( 1+a \ This 
issue has motivated the consideration of truncated Levy processes that exhibit long range 
dynamics while preserving the finiteness of some 24j or all moments 1_9|, |20|, |23[]. In the 
case of exponentially truncated processes, the Levy density decays as ~ x~^ l+0 ^e~ Xx , and 
for A 7^ all the moments are finite. In the present paper we limit attention to this type 
of processes and study the truncated fractional Fisher-Kolmogorov in which the Laplacian 
operator in Eq. (CQ) is replaced by the truncated fractional diffusion operator proposed in 
Ref. [23 1 . 

The organization of the rest of the paper is as follows. The next section reviews material 
on truncated Levy flights, defines the A-truncated fractional derivatives and discusses the 
fundamental solutions and scaling properties of the truncated fractional diffusion equation. 
Section III contains the core of the numerical results obtained from the integration of the 
truncated fractional Fisher-Kolmogorov equation for front-type initial conditions. Section IV 
presents an analytical asymptotic study based on the leading edge approximation. The 
conclusions are presented in Sec. V. 
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II. FRACTIONAL DIFFUSION WITH TRUNCATED LEVY FLIGHTS 



In the standard fractional diffusion model the transport of a scalar is governed by the 
equation 

c DU(x, t) = -ad x <j> + c [I ^D a x + r X D^] , (2) 

where we have included in addition to the spatial and temporal fractional operator an 
advective term. The operator on the left-hand side is the regularized (in the Caputo sense 



251 ]) fractional time derivative 



with < (3 < 1 and the operators on the right hand side of Eq.([2]) are the left and right 



Riemann-Liouville fractional derivatives 



25 



26] 



1 d ri 
r(m - a) dx m J a (x - y)<*+ l - m 



- D x<P=TT7- -T^T / 7Z „.wi- m d y» ( 4 ) 



(-l) m d m r b (f) 
tD ^ = T(m-a)dx^ L {y-xY+ 1 - mdy ' (5) 
with m — 1 < a < m. The weighting factors I and r are defined as 

i = — 2-^- r = _ + ( 6 ) 

2cos(a7r/2) ' 2cos(a7r/2) ' v ; 

with — 1 < < 1. According to Eqs. and (J6]), the parameter 6 determines the degree 
of asymmetry of the fractional operator. For 6 = the contributions of left and right 
derivatives are equal and the operator is symmetric. In the extremal, fully asymmetric case 
of main interest in this paper, 9 = — 1, and only the left derivative is present in the diffusion 
operator. Equation describes the fluid limit of a continuum time random walk (CTRW) 
in the case when the waiting time distribution function exhibits algebraic decay of the form, 
ip ~ and the particle jumps follow an a-stable Levy distribution, see for example [f| 

and references therein. 

The fractional equation (j2J) has found applications in several areas of physics, engineering 
and biology. For a recent discussion on the basic theory and applications of fractional 
diffusion see [5j. However more general transport equations can be obtained when a wider 
class of stochastic processes are considered. In particular, in the case when the particle 
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jump probability density function corresponds to an exponentially truncated (tempered) 



distribution with characteristic exponent of the form 



20] 



(1 + 0)( A + ik)« + (1 - 0)(A - ik) a - 2A«, 

Aet = ^ r7 < (7) 

2 cos(a7r/2) I (1 + 0)(A + ik) a + (1 - 0)(A - - 2A Q - 2ika9\ a - 1 , 

for < a < 1 and 1 < a < 2 respectively, the fluid limit of the corresponding CTRW leads 
to the equation [23j 

c dU = -Vd x <f> + cv:> x 4> - n4> , (8) 

where the A-truncated fractional derivative operator of order a, T>"' A , is defined as 

V a / = le~ Xx e Xx + re Xx X D^ e~ Xx . (9) 

The effective advection velocity is V = a for < a < 1, and V = a — v for 1 < a < 2 with 

" = icos(W2)r (10) 

and 

V= 7— r . (11) 

P cos(avr/2) v ' 



According to Eq. ( jTOI) . in the case 1 < a < 2, the truncation gives rise to a drift that 
depends on the asymmetry of the process. The parameter A determines the truncation of 



the tempered Levy process whose corresponding Levy density is given by 20j 



c I!+a e -A|*| forx<0, 

c^a;-( 1+a )e- Ax for 2 > 0, 

< a < 2, c> 0, -1 < 6> < 1 and A > 0. As expected, for A = 0, Eq. (EE]) reduces to the 
a-stable density, and Eq. ([8]) reduces to Eq. (J2]). 

The general solution of Eq. (jSJ) for an initial condition <p(x,t = 0) = <p(x, 0) is 

/oo 
G aj p,e,\(x - x', t)cj)(x', 0)dx' , (13) 
-oo 

where G a ,p,o,\ is the Green's function or propagator which corresponds to the solution with 
initial condition <f>(x,t = 0) = 5(x). The sub indices of G state explicitly that in general the 
solution depends on four parameters: the order of the fractional derivative in space a, the 
order of the fractional time derivative /3, the asymmetry of the fractional operator 9, and the 
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truncation A. Using the Fourier transform properties of the truncated fractional derivative 
it follows that 

1 r°° r i 

G aMX (x, t) = — e~ ikx E p \t^K ET {k)\ dk , (14) 

where Ep denotes the Mittag-Leffler function and Aet is given in Eq. (J7J). At first sight 
the linear term, — u<j), on the right hand side of Eq. ([8]) seems strange and likely to give 
rise to a unphysical damping of the transported field, 0. However, quite to the contrary, 
this term is critical to guarantee the conservation of (p. When <j) is interpreted as a prob- 
ability density function, this term guarantees the normalization and conservation of the 
total probability. One way to see this is to note that this term comes from the term pro- 
portional to — 2X a in (JTj) that implies AEr{k = 0) = 0. According to Eqs. ( 1T31) and (f!4l 
/ 4>(x,t)dx = Ep t^KErik = 0) J <po(x)dx which guarantees the conservation of (j> provide 
A ET (k = 0) = 0, since Ep(0) = 1. 

In this paper we focus on the special case of Eq. ([8]) with 1 < a < 2, (3 = 1, and 9 = — 1. 
Also, we assume an advection velocity a = v , which results in the following asymmetric, 
truncated fractional equation 

d t <j> = X [e- Xx -ooD a x (e A >) - A>] , (15) 

with x — c l |cos (ct7r/2)|. We restrict attention to this special case because it corresponds 
to the truncated version of the a-stable asymmetric fractional operator used in the front 
acceleration problem discussed in Ref. [h]]. For this case, the general solution in Eq. (fT4"|) 
reduces to 

G a ! _x A = ±- r e~ lkx+ ^ x -^ a - x ^ dk , (16) 

Z7T J-oo 

which can be equivalently written as 

GW-i-a = e- x *-* xat ( X ty 1/a G Q) i,-i,oM , (17) 

where 

1 r°° 

G a ,i,- lfi (7i) = — j^+^dk, (18) 



2lT J-oo 

is the Green's function of the asymmetric, a-stable (A = 0) fractional diffusion equation, in 
terms of the similarity variable 77 = x (xt)~ 1//a . From here, using the asymptotic expression, 



G a ,i-i,o(v) ~ V'^ a for 77 > [23, |2a] it follows that 



-Xx 

G^x-xte-*^ — , (19) 



for x > and x ^> For the decay of the left tail, we use the asymptotic expansion 

G^i -1,0(77) ~ \ri\ a2 e~ b2 ^ C2 for rj < and \r)\ > 1, where a 2 = (2 - a)/(2(a - 1)), b 2 = 



(a - lja"/^ 1 ) and c 2 = a/(a - 1) 28], to conclude 



G a ,i,-i,A ~ (Xty (a2+1)/a e~ xXat \x\°* e~ b ^~ C2/a , (20) 

for x < and |x| ^> (xt) 1 ^". Since we are assuming that 1 < a < 2, it follows that 
C2 > 1 and the — \x\ C2 term in the exponent dominates the X\x\ term, leading to a faster 
than exponential decay of the left tail for any value of A. Figure 1 shows plots of the 
Green's function in Eq. (|16|) for a = 1.5 and different values of A, along with the asymptotic 
approximation in Eq. f|T9|) . 

An important property of truncated Levy flights, originally discussed in Refs. 0,0, EL 



is the ultraslow convergence to Gaussian statistics. According to this result, the crossover 
time for Gaussian behavior to appear, r c , scales as 

T C ~ X - 1 \~\ (21) 

as expected, as A — > 0, r — > 00. When memory effects are incorporated using fractional time 
derivatives, the crossover dynamics is richer. In particular, for 2/J/ck > 1, r c ~ %" 



231 ] . The time scale r c 



signals the crossover from super diffusive to subdiffusive dynamics 
will play an important role in the dynamics of the fronts. 



III. FRONT PROPAGATION IN THE PRESENCE OF TRUNCATED LEVY 
FLIGHTS: NUMERICAL RESULTS 

To study the role of truncation in the superdiffusive acceleration of fronts due to 
Levy flights we consider the fractional Fisher-Kolmogorov equation originally introduced 
in Ref. (lOj and substitute the Riemann-Liouville fractional derivatives (which correspond 
to a-stable Levy processes) by the truncated fractional derivative in Eq. In the most 
general case the resulting equation is 

C DU = -Vd x <f) + cV a /<p - ^ + 70 (1 - 0) . (22) 



However, as mentioned before, to compare the results with those in Ref. 10J we will restrict 



attention to asymmetric, truncated fractional diffusion operators of the form in Eq. ffT5l) and 
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consider 

d t cf> = X [e~ Xx -ooD: (e Xx 4>) ~ A>] + T (1 - 0) . (23) 

In this section we present results obtained from the numerical integration of Eq. fl23|) . We 
assume = A, for x < 0, where A is a constant and discretize the fractional derivative in the 
x G (0, 1) domain using the Grunwald-Letnikov representation. Details of finite-difference 



methods for the solution of fractional diffusion equations can be found in |29l . |30j. In all 
the numerical simulations we consider a = 1.5, 9 = —1, 7 = 1, x = 5 x 1CT 7 , and initial 
conditions of the form 

= 0) = e" ra , (24) 

where v is a constant. 

Figure [2]- (a) shows snapshots of the front profile as function of x at different times in 
the a-stable (A = 0) case obtained from the numerical solution of Eq. ff23l . In this case the 
: ront exhibits an algebraically decaying tail which in log-log scale manifests as a straight line 
lo| . Figure 12]- (b) shows that the algebraic decay of the tail remains for small values of A. 
In fact, as we will discuss in the next section, there is an intermediate asymptotic regime in 
which the role of truncation is negligible. Outside the intermediate asymptotic regime the 
effect of the truncation depends critically on the ratio of the length scale of the truncation, 
1/A, and length scale, 1/V, of the initial condition. When, A < v, i.e., when the initial 
condition decays faster than the truncation, the tail of the front scales as ~ x~ x ~ a e~ Xx 
as shown in Fig. G9-(c). On the other hand, when A > u, the truncation effects dominate 
and Levy statistics seems to have no effect on the dynamics of the front which, as shown in 
Fig. [2]-(d) (note the log- normal scale), exhibits the usual exponential decay of the diffusive 
Fisher-Kolmogorov model. 

When the front exhibits "rigid" propagation with a constant velocity, as in the diffusive 
Fisher-Kolmogorov case, it is straightforward to define and numerically compute the front 
speed. However, when the front accelerates and deforms, as it is the case in the a-stable 
and truncated fractional Fisher-Kolmogorov equation, computing the front speed is not 
straightforward. One way to approach this problem is to consider the Lagrangian trajectory 
of the front's tail. For a given value of <fro, we define the Lagrangian trajectory, xl = xi{t), 
of the front according to cj)(xi,(t),t) = 4>q. Given x L we define the Lagrangian velocity as 
vl = dxi/dt and the Lagrangian acceleration as = dvi/dt. 
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Figure [3] shows space-time plots of the Lagrangian trajectories of fronts corresponding 
to 0o = 10~ 6 . Due to the relatively small value of <po, these orbits follow the Lagrangian 
dynamics of the fronts' leading edge. The solid lines denote the numerical values for different 
values of A and the dashed lines denote the result of the asymptotic analytic calculation that 
will be discussed in the next section. For < A < 1 it is observed that the front propagates 
very fast. However, as A increases the speed of the front is reduced. The corresponding 
Lagrangian velocities and accelerations are shown in Figs. H] and [51 For small A, the font 
acceleration grows monotonically. However, for larger values of A the Lagrangian velocity 
approaches a terminal velocity as t — > oo. Note that as A increases, the terminal velocity is 
smaller and the convergence is faster. 

As shown in Fig. [5J in all cases the acceleration exhibits a pulse like behavior in the time 
evolution. The time of peaking of the pulse increases (approximately exponentially) with A 
and the amplitude of the acceleration's peak decreases (approximately exponentially) with A. 
However, the key feature to note is the evolution of the acceleration following the transient 
short pulse. For small A the front acceleration exhibits a monotonic increase whereas for 
larger values the acceleration exhibits a transient growth followed by an eventual decay. 
Figure [6] shows the asymptotic scaling behavior of the Lagrangian velocity and acceleration. 
The solid line curves denote the numerical results for different values of A. It is observed 
that the convergence to the terminal velocity scales as v* — Vl ~ 1/t and the decay of the 
transient acceleration scales as ol ~ 1/t 2 . These numerical results are in agreement with the 
asymptotic scaling, shown with dashed line curves, that will be discussed in the following 
section. 

IV. LEADING EDGE ASYMPTOTIC DYNAMICS: ANALYTICAL RESULTS 

In this section we present analytical results describing the asymptotic behavior of fronts in 
the truncated fractional Fisher-Kolmogorov equation. The results are based on the leading 
edge approximation. This approximation exploits the idea that at the leading edge of the 
front <fi <C 1, and therefore in this region the nonlinear reaction term can be linearized 
around = 0, resulting in the linear equation 




(25) 
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Note that the truncated fractional derivative has a direct effect on the growth rate through 
the term — x^ a - As discussed in the previous section, this term is key to guarantee the 
conservation of the transported field. Two characteristic time scales can be distinguished 
in this problem: the cross-over to Gaussian statistics time scale, r c ~ 1/ (xX a ), and the 
reaction time scale, r r = I/7. We will assume that r c > r r to guarantee that the effective 
reaction constant, j e ff = 7 — x^", is positive as needed for the excitation and propagation 
of "pull" type fronts. 
Substituting 

= e- Ax+{7 -* AQ) V(a;,t), (26) 
into Eq. (|25p gives the asymmetric, a-stable fractional diffusion equation 

d^ = X -ooD^. (27) 

The general solution of this equation for an initial condition ip(x,t = 0) = i(jq(x) can be 
written as 

f°° * T l/a 1 

ip(x,t) = / ^,1-1,0(^)^0 x-(xt) /a rj drj, (28) 
J— 00 L J 

where G a) 1,-1,0 given in Eq. (|T8l) . 

Consistent with the numerical simulations, we consider an initial condition of the form 
<ft(x, t = 0) = A for x < and <fr(x, t = 0) = e~ ux , where A and v are constants. Substituting 
the corresponding initial condition for 1/1, according to Eq.( l26l) . into the solution in Eq.(l28l) 

we get 

rx/r A poo A 

i> = e-^~ x)x / e^-^G^Mdri + Ae Xx / G Q , 1 ,_ 1 , (r ? ) e - A ^ , (29) 

J — 00 Jx/t 

where we have defined r = (xt) 1 ^ . In terms of <p the solution can be written as 

= e —+h-x^)t Ix + Ae^- xXa H 2 , (30) 

where 



/ x l r ~ f°° ~ 

e^- x ^G aA ^ ( V )d V , 1 2 = / G ajl) _ 1)0 (7 7 ) e - AT ^7 7 . (31) 
-00 Jx/t 



r 1 — "a,.,.-..o(//)'///. J2 ■ , 

/x/r 

The goal of this section is to study the asymptotic behavior of T\ and X 2 for x/r — > 00. 
Before getting in the calculation of main interest here, it is instructive to consider first the 
diffusive (Gaussian) and fractional (a-stable) limits of the leading edge solution in Eq.fl30l). 
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A. Diffusive case 



In the diffusive (Gaussian) case, (a, 6, A) = (2, 0, 0), the leading edge solution in Eq. (J30 
reduces to 



/ X /Vxt * , roc 

G 2AA0 ( V y^t dv + Ae yt / G 2 , 1A0 dr, 
-oo Jz/Vxi 



(32) 



where 62,1,0,0 — Vl^V^") 6 ^ ^ i s the Gaussian propagator. Introducing the normal proba- 
bility distribution function P(z) = {l/y/Zn) Jl^e^^^du, 



P 



. ) +Ae 



i-p 



X 



(33) 



Using the asymptotic expansion P{z) ~ 1 — (l/v / 27r)e z2 ^ 2 /z, we conclude that in the limit 
77 = x/t 1, x > 2uxt, 



71 



A 



where 



x x — Ivyt 



T 1 

c = - + v\ ■ 
v 



-(x-c m i)(x+c m t) 



(34) 



(35) 



That is, in this case, the leading edge exhibits the well-known asymptotic exponential de- 
pendence, <p ~ e - u ( x ~ ct )^ anc l the front propagates at the constant speed c with c m = 2^Ffx 



corresponding to the minimum speed achieved for v = J^f/x- Note that according to 
Eq. (1341) in this case the convergence to constant speed is exponentially fast, i.e., the second 
term in the asymptotic expansion scales as ~ e~ a { x / T ) 2 . This result will be contrasted below 
with the much slower convergence in the case of truncated Levy flights. 



B. Fractional case 

The fractional (a-stable) case was discussed in Ref. |lo| . In this case, the leading edge 
solution is given by Eqs. (l3"0l) and ( 13~TT) with A = and the corresponding leading asymptotic 
behavior is 



Xte 



7* 



A 1 

-x~ a + -x~ l ~ a + . 



(36) 



.a v 

where as mentioned before the constant A relates to the boundary condition <p = A for 
x < 0. The critical difference with the Gaussian case is the algebraic decay of the leading 



12 



edge accompanied by the exponential acceleration of the front. Note that when A ^ 0, 
(which was the case considered in Ref. [h]]) the front tail exhibits the decay <\> ~ l/x a . 
However, when A = 0, the front decays faster, ~ l/x a+l . Figure [2]-(a) shows a numerical 
verification of this scaling. This result will be contrasted with the truncate Levy flights case, 
where as in the Gaussian case, it will be shown that the rate of decay of the front's tail is 
independent of A. 

In the Gaussian case, the spatio temporal evolution of the leading edge depends to leading 
order on the variable x — ct which implies a "rigid" translation of the exponential tail of 
the front and allows the interpretation of c as the front. However, in the non-Gaussian case 
each point of the leading edge moves at a different speed and the tail does not translate 
rigidly. As discussed in the previous section we circumvent this problem by considering the 
Lagrangian trajectory xl = xi(t; O ) of a point in the leading edge of the front such that 
0(xi(t),t) = 0o where O 1. According to Eq. (|36|) . in the a-stable case, for A = 0, 

1 



xiit) = C exp 



(7* + lnt) 



.1 + a 

where C is a constant that depends on O , \ an d v, and 



(37) 



VL{t) 



1 



eW=)*, (38) 



1 + a \0o^ / 

which implies an unbounded, exponential growth of the front speed. For large t, the corre- 
sponding leading-order behavior of the front acceleration is 

^(W^f^l^eteK (39 ) 



1 + aJ V^o 2 ^ 

As shown in Figs. [3HU and [5] the asymptotic results in Eqs. fl37|) - fl39l) are in good agreement 
with the numerical results discussed in the previous section. 



C. Truncated case 



Going back to the general truncated fractional case, we consider first the asymptotic 
behavior of I2 in Eq. (13TT) . In the limit x/t — » 00 , the integration variable satisfies 77 ^> 1 
and thus we can use the asymptotic expression of the Green's function corresponding to the 
right, algebraic decaying tail, G a ,i,-i,o ~ and get, after an integration by parts ,the 
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asymptotic expansion 

/•OO 

X2 ~ / 1 

Jx/t 



drj 



A x 1+a 



a + 1 
At A/r 



Arr?^-(a+2) 



Integrating by parts once more, gives the next term in the asymptotic series 



1-2 



t e 



-Xx 



A x 1+a 



(a + 1) 
Ax 



+ . . . 



To deal with X 1; introduce a cut-off Q such that 1 <C Q and write the integral as 

fx/ 



li = C+ [ e^^G^Mdrj. 



(40) 



(41) 



(42) 



where the constant on the right hand side is defined as C = e( u ~ x ^ TV G at i^ifi(r])dr]. 
Note that because of the faster than exponential decay of the asymmetric Levy distribution 
Ga.i-i.ofe) f° r V < in Eq. pO]) the integrals converge for either sign of v — A. Since both 
the cut-off Q and x/r are assumed to be large, we can substitute as before the asymptotic 
expression of the Green's function in the integral of Eq. fl42p . and after an integration by 
parts obtain the asymptotic expansion 



^(y— \)x 



+ . 



v - A x 1+a 

Substituting Eq.g3]) and Eq.flilj) into Eq.(pl) we get 

\u-X XJ x a+1 



(43) 



(44) 



for v 7^ A. The issue now is to determine the leading order term in Eq. (jUJ). The answer to 
this problem depends on the relative values of v and A. 

If v > A > 0, i.e., if the initial condition decays faster than the truncation, the leading 
order term in Eq. (|44p for large x is 

tx 



1 A 
u-X + J 



-\x+{~f~ X ^ a )t 



-1+a 



(45) 



Note that, contrary to the a-stable case, the asymptotic spatial decay of the front leading 
edge in Eq. (145!) is independent of the value of A. The role of the truncation is clearly seen 
in the exponential factor e~ Xx that dominates the decay for x ^> 1/A. Figure 2-(c) shows a 
very good agreement between the numerical result and the scaling in Eq. ( ]4"5l) . 
When 1 < 1/A we can expand the exponential in Eq. (|45p and write 



1 A 
v- X + T 



txe 



(7-xA a )« 



x 



l+a 



X A ' 2 

Xx + —x 



(46) 
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According Eq. (|46[) in the intermediate asymptotic regime, (Xt) 1/a < x < 1/A, the front 
exhibits to leading order the ideal (untruncated) Levy flight algebraic scaling, ~ l/x 1+a . 
This scaling is numerically verified in Fig. 2-(b). Moreover, as Figs. 3, 4 and 5 show, in the 
intermediate asymptotic regime (which in the numerical simulations roughly corresponds to 
< A < 1) the front's velocity and acceleration exhibit unbounded monotonic growth and 
follow to a good approximation the ideal Levy flight scaling in Eqs. f l3"T]) - ([3"9"]) . 

Going back to Eq. (j4"5]) we have that outside the intermediate asymptotic regime, i.e. for 
x > 1/A, the Lagrangian trajectory of the front is given by 

-Xx L {t) + (7 - xA a )t + hit - (a + l)\nx L (t) = M , (47) 

where M is a constant that depends on 0o- From Eq. (1471) we obtain the following expression 
for the Lagrangian velocity of the front, vl = dxL{t)/dt, 

Mt) = 7 * 7* • (48) 

A + 

As shown in Figs. [3] and H] there is good agreement between the numerical results (solid 
lines) and Eqs. ( HTl) and ( 1481) (dashed lines) in the asymptotic regime x ^> (xt) 1 ^". From 
Eq. (j48p it follows that in the limit t — > 00 

^ L ~ W * ~ At " ' ' ^ 



where the terminal velocity is given by 



(50) 



A ' 

which is positive since it has been assumed that 7 > X a x- The asymptotic approach to the 
terminal velocity is clearly observed in Fig. H] where the horizontal dashed lines show the 
terminal velocity in Eq. (1501) for the values of A considered in the numerical simulations. 
According to Eq. (149]) the time required for the front velocity to approach the terminal 
velocity within a given margin v* — vl scales as t ~ 1/A. The corresponding Lagrangian 
acceleration of the front, = dvi(t)/dt, is given by 

1 2 



adt)= {in 1) 

t (Ati>* + 1) 



v L (t) 



x L (t)/t 

From Eqs. (148]) and (I5T]) it follows that for large times, 



(51) 



«L(t)~^. (52) 
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As shown in Fig. [6j the analytical scaling relations agree well with the numerical results. In 
this figure, the curved dashed lines correspond to the analytical result in Eqs. (J15]l and fl5"Tj) 
and the straight dashed lines correspond to the scaling in Eqs. (15"Uj) and (I5"2j) . Thus, outside 
the intermediate asymptotic regime, i.e. for 1/A < x, the front acceleration decays and the 
front approaches a constant terminal velocity as t — > oo. However, the convergence of the 
dynamics to the constant front velocity regimen exhibits a very slow, ~ l/(At), algebraic 
decay compared to the significantly faster exponential convergence in the diffusive case. 

The calculations presented up to now assumed v > A. However, when A > u, i.e., when 
the initial condition decays slower than the truncation, the leading order term in Eq. (jHJ) is 

~ Ce-^-^ , (53) 

which indicates that the front moves with the constant velocity 

v=l^p<. (54) 

An example of this over-truncated regime is presented in Fig. |2]-(d) that shows a very good 
agreement with the exponential decay in Eq. (1531) for A = 100 and v = 50. The remaining 
case to consider is A = v. In this case, the asymptotic approximation in Eq. (|41[) still holds. 
However, the expression in Eq. (TJ3]) can not be used, and we have to go back to the integral 
Xi Eq. flUD 

Xi = C + £ /T ^,1,-1,0^)^, (55) 

where this time the constant on the right hand side is defined as C = G a ,i~i,o(v)df]- 
Using the asymptotic expression for G ajl _ li0 ~ rf x ~ a and integrating, it is concluded that 
the leading order term in Eq. fl3"U|) is (ft ~ c e - Ax '+(7-x^ Q: )* w hich implies that the constant 
velocity of the front is given by Eq. fl54|) with v = A. 



V. CONCLUSIONS 



We have presented a numerical and analytical study of the role of truncated Levy flights 
in the propagation of fronts in reaction superdiffusion systems. The study was based on the 
truncated fractional Fisher-Kolmogorov model in which the spatial derivative is replaced by 
the A-truncated fractional derivative of order a. 
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Depending on the level of truncation four front propagation regimes can be distinguished: 
an asymptotic algebraic regime, an intermediate asymptotic algebraic regime, a truncated 
regime, and an over-truncated regime. The asymptotic algebraic regime corresponds to 
A = 0. In this case the problem reduces to the a-stable fractional Fisher-Kolmogorov 
problem that exhibits exponential acceleration and algebraic decaying tails for x ^> (xt) 1 ^". 
The intermediate asymptotic regime is characterized by (xt) 1 ^* <C % <C 1/A. We have 
shown numerically and analytically that in this regime the truncation effects are negligible 
and the algebraic decay of the tail as well as the acceleration of the front prevail. Outside 
the intermediate asymptotic regime, i.e. for x > 1/X and (xt) 1 ^ "C x, the tail of the 
front exhibits the tempered decay <fi ~ e~ Xx /x 1+a , the acceleration is transient, and the 
front eventually reaches a terminal speed as t —>■ oo. In the over-truncated regime the 
truncation decays faster than the initial condition, i.e. A > v. In this case, Levy flights 
have apparently no qualitative effect on the asymptotic dynamics of the front that exhibits 
a diffusive-type exponential tail, ~ e~ ux , and constant propagation velocity, v = j e ff/is. 
However, contrary to the diffusive case, the constant velocity in the over-truncated case is 
a monotonically decaying function of v and has no finite minimum. 

Although in the truncated regime the acceleration decays and the front eventually reaches 
a constant terminal speed, the numerical and analytical results show that the convergence 
of this process is very slow. In particular the front acceleration asymptotically decays as 
ol ~ ot/ (At 2 ) and the approach to the terminal velocity scales as vl ~ f* — a/ (At). This 
algebraic convergence is in sharp contrast with the exponential convergence observed in 
the diffusive Fisher-Kolmogorov equation. In this sense the truncated regime resembles the 
"ultra-slow" convergence regimen observed in transport problems without reaction terms. 

One of the motivations and potential applications of the present work lays on the study 
of transport in magnetically confined plasmas. In this system it has been suggested that 
reaction diffusion models provide insightful, though highly simplified, reduced descriptions 
of the interaction of turbulence and transport. Current models usually assume Laplacian 
diffusive operators, see for example Refs. |31|. However there is evidence that transport in 
magnetically confined plasmas is not necessarily diffusive, see for example and references 
therein. On the other hand, it has been argued that truncated Levy distributions might play 



a role in the description of electrostatic potential fluctuations 



331 ] . Based on this, it would 



of interest to explore the implications of the present work on the corresponding plasma 
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transport models. 

Throughout this paper we have limited attention to extremal, 6 = — 1, transport processes 
and assumed an external advection velocity to cancel the drift resulting from the truncation 
of the asymmetric fractional derivative. A problem of interest is to perform more general 
numerical simulations considering different degrees of asymmetry and including memory 
effects through the use of fractional derivatives in time. As mentioned in the introduction in 
recent years several papers have discussed the role of Leevy flights in front propagation and 
pattern formation in reaction-anomalous-diffusion systems. It would be of interest to explore 
the role of truncation effects in these systems. Beyond its intrinsic theoretical interest, the 
study presented here might have relevance in the interpretation and modeling of laboratory 
experiments and numerical simulations of complex systems. It is plausible that in these 
systems, the presence of boundary conditions, finite size domain and decorrelation could 
introduce "non-ideal" Levy-flights dynamics of the type discussed here. 
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FIGURE CAPTIONS 



FIG. 1. (Color online) Green's functions of the asymmetric truncated fractional diffusion 
Eq. ffTBT) for a = 1.5, 9 = —1, and A = 0, 10, and 20. The solid (blue) lines show 
G as function of x for fixed t, and the dashed (black) lines show the corresponding 
asymptotic dependence according to Eq. (fl9|) . The numerical and asymptotic re- 
sults are practically indistinguishable for x > 0.1 and for visualization purposes the 
asymptotic result has been shifted downward a little bit in the plot. 

FIG. 2. (Color online) Depending on the value of the truncation parameter A four front prop- 
agation regimes can be distinguished in the truncated fractional Fisher-Kolmogorov 
Eq. fl23l) : (a) Asymptotic algebraic regime for A = 0; (b) Intermediate asymptotic 
algebraic regime for A 7^ 0; (c) Truncated regime for < A < v\ (d) Over-truncated 
regime for A > v. In all four panels the left most (red) curve denotes the initial 
condition in Eq. (|2~4"|) and the other (blue) solid line curves show the front profiles 
at different times. The dashed (red) curve shows the analytical asymptotic result 
according to Eq. (1361) with A = in panels (a) and (b); according to Eq. (j4"5"l) in panel 
(c); and according to Eq. (1531) in panel (d). In all cases a = 1.5, 9 = —1, 7 = 1, and 
x = 5 x 10~ 7 . 

FIG. 3. (Color online) Space-time Lagrangian front paths, 0(xi(t),t) = 0o with O = 10~ 6 , 
according to the asymmetric truncated fractional Fisher-Kolmogorov Eq. ff23l with 
a — 1.5, 9 — —1, 7 = 1, x — 5 x 10~ 7 and different values of A. The solid (blue) 
curves denote the numerical results and the dashed (red) curves the asymptotic result 
according to Eq. (H7|) . The dotted (green) lines denote the front speed (upper line) 
and the minimum front speed (lower line) in the Gaussian diffusive case according 
to Eq. fl35l) . As the A = 1 case shows, the effect of truncation is negligible in the 
intermediate asymptotic regime. 

FIG. 4. (Color online) Time dependence of Lagrangian front velocities, vi(t) = dxL(t)/dt, 
according to the asymmetric truncated fractional Fisher-Kolmogorov Eq. (!23j) with 
a = 1.5, 9 = —1, 7 = 1, x = 5 x 10~ 7 and different values of A. The solid (blue) 
curves denote the numerical results and the dashed (red) curves the asymptotic result 
according to Eq. (HBl) . The horizontal (black) dash lines denote the corresponding 



21 



terminal velocities according to Eq. (1501) . The solid (green) line at the bottom 
denotes the corresponding front speed in the Gaussian diffusive case according to 
Eq. dSSD- 

FIG. 5. (Color online) Dependence of front acceleration on truncation. The panels show the 
time dependence of the Lagrangian front acceleration a^(t) = dvL(t)/dt, according 
to the asymmetric truncated fractional Fisher-Kolmogorov Eq. ( 1231) with a = 1.5, 
9 = — 1, 7 = 1, x = 5 x 1CT 7 and different values of A. The solid (blue) curves denote 
the numerical results and the dashed (red) curves the asymptotic result according 
to Eqs. (1391) and fl5T]) . In the intermediate asymptotic regime, panels (a)-(c), the 
front exhibit monotonically increasing acceleration. In the truncated regime, panels 
(d)-(f) the acceleration is transient and decays at large times. 

FIG. 6. (Color online) Asymptotic algebraic scaling of approach to terminal velocity, v* — vl, 
and asymptotic algebraic scaling of Lagrangian front acceleration decay . The solid 
(blue) curves denote the numerical results for a = 1.5, 9 = —1, 7 = 1, x — 5 x 10~ 7 
and different values of A. The dashed (red) curves show the asymptotic results 
according to Eqs. (j48p and (1511) . and the dashed (black) straight lines the analytical 
asymptotic scalings, v* — vl ~ a/(Xt) and ~ a/(\t 2 ) according to Eqs. (Pf9l) and 

(EH). 
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